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I. INTRODUCTION 

An important point in the dynamics of open quantum 
systems is the initial state preparation. Due to unavoid¬ 
able interactions, there generally exist initial correlations 
between a system and its environment. Thus any physi¬ 
cal process of preparation of the initial state of the system 
will affect the state of the environment as well. The ques¬ 
tion then arises of how the preparation procedure and 
initial correlations influence the subsequent evolution of 
the open system. Different aspects of this question were 
discussed by many authors (see, e.g., Refs. [ll-Q and ref¬ 
erences therein). Of special interest is the decoherence 
phenomenon (the environmentally induced destruction of 
quantum coherence). For instance, decoherence plays a 
crucial role in the dynamics of two-state systems (qubits) 
which are the elementary carries of quantum informa¬ 
tion 

At first sight it is natural to expect that initial cor¬ 
relations between an open quantum system and an en¬ 
vironment with a huge (or even infinite) number of de¬ 
grees of freedom would increase the decoherence rate. 
This is indeed the case for some preparation proce¬ 
dures iii but is not true in general. We refer to the 
paper , where the authors demonstrated with examples 
of specific qubit-environment models that for some ini¬ 
tial system-environment correlated states, the “purity” 
of the qubit is greater than in the case of initially uncor¬ 
related states. Although the results of Ref. Q rely on a 
somewhat artificial assumption that the environment is 
initially prepared in a pure quantum state, the fact that 
quantum coherence in open systems can be enhanced due 
to system-environment correlations seems as itself to de¬ 
serve thorough studies. It could provide, for instance, 
a way for constructing dynamics of qubit registers with 
interesting and very promising properties. 

In this paper we aim to examine in detail the con¬ 
nection between the initial preparation and the qubit 
dynamics. In particular, this allows us to establish the 
conditions for appreciable environment-induced enhance¬ 
ment of quantum coherence in a qubit. A somewhat un¬ 
expected result is that the purity of a qubit state can even 


increase with time. The principal difference between our 
analysis and that presented in Ref. @ lies in the inter¬ 
pretation of initial states of the composite (qubit plus 
environment) system. As we have already mentioned, 
in Ref. @ the initial states were taken in the form which 
could illustrate the role of qubit-environment correlations 
but, unfortunately, these states are very unlikely to be 
physically realizable. In the present paper we consider 
more realistic preparation procedures based on quantum 
measurements pI]Lflll|. 

The structure of the paper is as follows. In Sec. m 
we give a brief review of measurement schemes for open 
quantum systems and applications to qubits. In Sec. IIIII 
we treat time evolution of a qubit coupled to a bosonic 
environment through a dephasing interaction. Exact ex¬ 
pressions are given for the time-dependent elements of 
the qubit density matrix when the initial state is pre¬ 
pared by a selective or non-selective measurement. Our 
central goal in this section is to demonstrate that for 
a large class of initial states prepared by non-selective 
measurements, the coherences (off-diagonal elements of 
the qubit density matrix) increase with time at the ini¬ 
tial stage of evolution. In Sec. IlYl we consider the time 
behavior of the qubit purity and entropy in the regime 
of enhancement of coherence. Conclusions are drawn in 
Sec. El 


II. SELECTIVE AND NON-SELECTIVE 
PREPARATION MEASUREMENTS 

As an introduction to our subsequent development, we 
start with a brief discussion of a rather general quan¬ 
tum measurement scheme which can be used to construct 
statistical ensembles describing initial states of real open 
systems. We then apply this scheme to a qubit interact¬ 
ing with its environment. 

Suppose that at all times 0 < t an open system S is 
in thermal equilibrium with a heat bath B, and at time 
zero one makes a measurement on the system S only. 
AccordiM to g eneral principles of quantum measurement 
theory |10l-[l^. the state of the composite system {S + B) 
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after the measurement is described by the density matrix 
Sssi^) = ^mgeg^Li (1) 

m 

where Pgq is the equilibrium density matrix at temper¬ 
ature T. Operators act in the Hilbert space of the 
system S and correspond to possible outcomes m of the 
measurement. In a particular case of a selective measure¬ 
ment, the system S is prepared in some pure state l"^). 
Then the sum in Eq. o collapses into a single term, so 
that 

(^) 

where is the projector onto the quantum 

state l"^) and Z is the normalization factor. In general, 
the density matrix © describes the resulting ensemble 
after a non-selective measurement in which the outcome 
m may be viewed as a classical random number with the 
probability distribution 

w{m) =Ti{F,^g^^} . ( 3 ) 

Positive operators are called the “effects”. 

Here and in what follows, Tr denotes the trace over the 
Hilbert space of the composite {S -I- B) system, while the 
symbols Tr_g and Tr^ will be used to denote the partial 
traces over the Hilbert spaces of the system S and the 
heat bath, respectively. In order that w(rn) be normal¬ 
ized to 1, the effects F^ must satisfy the normalization 
condition (the resolution of the identity) 


A few remarks are needed here. Whereas the states 1 ^/’^) 
form an orthonormal basis, for the transformed states 
\^Pra) this is true only if the unitary operators are 
identical for all outcomes m (i.e., = U). In such 

cases, since have, in addition to 

Eq. ([ 3 ]), another resolution of the identity 

( 9 ) 

m 

In other words, to the measurement scheme defined in 
terms of the effects ([6]) and the H-operators 0 there 
corresponds the “dual scheme” characterized by 

F^ = UF^U^ = \^^){ipJ, 

■ 

If the form of depends on the outcome to, then the 
transformed states |<^^) are not orthogonal in general. 
Some of these states may even be identical. 

Let us now apply the above general construction to a 
qubit. In the formal “spin” representation, the canonical 
orthonormal basis states of a qubit are 

!“>=(?)■ = ( 11 ) 

It is well known that all pure states |'i/'(a)) = |o) cor¬ 
respond to points of the sphere |d| = 1, where a = 
(0^,02,03) G Normalized state vectors are given by 
(see, e.g., Ol) 


Y,F^ = J 2 nln^ = i, ( 4 ) 

m m 

where / is the unit operator. 

The precise form of is determined by the details 
of the measuring device. We restrict ourselves to phys¬ 
ical situations in which some observable A with a dis¬ 
crete, non-degenerate spectrum is measured. Then 
we have [T 3 , [ill 

( 5 ) 

1 /2 

where Fm is the square root of F^, and the unitary 
operator describes the disturbance of the system S 
by the measurement device. Formula ([ 5 ]) is applicable 
even to approximate measurements where the spectrum 
is measured with finite resolution US nil. We will 
be content, however, with performing the analysis for the 
case of infinite resolution when the effects are written as 

Fm = l^m)(V'ml (6) 

and Eq. 0 yields 

( 7 ) 


with 


/ e cos(0a/2) \ 

V sin(0j2) y’ 


( 12 ) 


where (j)a and 9 a are the Euler angles of the unit vector 
a describing the direction of the “spin”. They satisfy 
Ol ia2 = sin 0 a = cos 0 a- The Euler angles 

corresponding to the state | — a) are 

0-a=T^- <l^-a = (l>a + '^- ( 13 ) 

Using these relations together with Eq. m gives 

/ sin(0a/2) \ 

\ cos(0a/2) J 


Linear operators in the qubit’s Hilbert space can be 
represented as linear combinations of the unity operator 
and the Pauli matrices. For example, the operator 


cr(d) = -I- (7202 + cr^a^ 


( 15 ) 


describes the spin component in the direction a. The 
state vectors | ± a) correspond to the eigenvalues ±1 of 
cr(a) and form an orthonormal basis, i.e., 


( 8 ) 


|d)(d| -k I -d){-d| = I. 


|(0m) = Uml'fpm) 


( 16 ) 
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Note that the states |a) and | — a) are related to the a qubit, which is described by Eqs. (l6])-(|8l) with different 
canonical basis states (HIJ by unitary operators Um- Recalling Eqs. m , we write 


\a) = U{a)\l), \-a) = U{dm ( 17 ) 

with the unitary operator 

/ g-i0a/2 cos(6»a/2) sin(6»a/2) \ 

U{a) = 

V sin(6l„/2) cos{9a2) J 

( 18 ) 

Selective measurements ([ 2 ]) on a qubit can be charac¬ 
terized by the projectors P{a) = |d)(d| while the general 
non-selective measurement scheme (0)-® is associated 
with three states: |d), |6i), and 162)- The effects and the 
ri-operators are defined as 

Fi = |a)(a|, ^^2 = I - a)(-d|, 

^i = \bi){a\, ^2 = \b 2 ){-a\- 

According to Eqs. (IT71) . we have 

\b,)=U{b„a)\d), \b^) = U{-b^,a)\-a), ( 20 ) 

where the unitary operator U ( 5 , a) is expressed in terms 
of the operators C®: 

U{h,d) = U{b)U\a). ( 21 ) 

Let us briefly consider some important special cases of 
the measurement scheme (HU): 

i) b^ = a, 62 ~ This is the simplest scheme cor¬ 
responding to = / in the general formula ®. Phys¬ 
ically, here we are dealing with non-selective measure¬ 
ments where the measuring device does not disturb the 
basis states |a) and | — a). In this case the fl-operators 
coincide with the effects: 


Ill = Fi = |a)(d|, O2 =F2 = I-a)(-d|. (22) 


Clearly, the same operators correspond to the dual mea¬ 
surement scheme m- 

ii) bi = b, b2 = —b, where b is an arbitrary unit vec¬ 
tor. This case corresponds to = U = U{b,a) in the 
formula ®. It follows from Eqs. (fl®) that 

El = |d){d|, F2 = |-d){-d|, 

( 23 ) 

IIi = |6)(d|, Il2 = I - ■ 

There exists the dual measurement scheme cni) with 


Fi = |6)(6|, F2 = |-6){-6|, 


hi) bi = b, b2 = b with an arbitrary unit vector b. This 
is an example of a non-selective measurement scheme for 


Ei = |a)(a|, F2 = |-d)(-d|, 

^i = l^)(a|> ^2 = \b){-a\. 

It is easily verified that 

|6) = ; 7 i|d), \b) = U 2 \-ci), ( 26 ) 

where Ui ^ U2 and are given by 

Ui = U{b, d), U2 = U{b, -d). ( 27 ) 

Note that in this case there is no dual measurement 
scheme. 

We close this section with a remark about the above- 
discussed non-selective measurement schemes. Math¬ 
ematically, all the schemes are generated by orthogo¬ 
nal resolutions of the identity o where the effects 
are projectors (F^ = F^) and are given by Eq. ®. 
This is a natural generalization of the well-known von 
Neumann-Liiders projection postulate for ideal quantum 
measurements (see, e.g., Ref. [la)- One can, however, 
construct more general measurement schemes associated 
with the notion of the positive operator-valued measure 
(POVM) [Hill. A POVM is defined by N positive op¬ 
erators F^ which form the resolution of the identity, but 
in general F^ ^ F^. Usually the operators F^ can be 
represented in the form ® where the set {|'0m)} over¬ 
complete, i.e., the number of outcomes N exceeds the 
rank of the density matrix of the system (for a qubit 
N > 2 ). In the present paper we will not consider 
such general situations and restrict ourselves to the von 
Neumann-Liiders projection measurements. 


III. THE DEPHASING MODEL: DYNAMICS OF 
DECOHERENCE 

A. Exact expressions for the coherences 

Our central goal in this section is to demonstrate that 
the qubit dynamics with initial states prepared by non- 
selective measurements exhibits a number of physically 
interesting and even somewhat unexpected features when 
compared with the case of selective measurements. This 
is especially important when one is dealing with decoher¬ 
ence phenomena. 

Our discussion will be based on the analysis of the 
simple dephasing model describing the main decoherence 
mechanism for certain types of system-environment in¬ 
teractions In this model, a two-state system 

(qubit) (S) is coupled to a bath {B) of harmonic os¬ 
cillators. Using the “spin” representation for the qubit 
with the basis states (ED, the total Hamiltonian in the 
Schrodinger picture is taken to be (in our units h = 1 ) 

H = Hs+Hb+ 

= yCTs + XI ^kblbk + 0-3 X( 5 fc^fc + 9 *kbk), ( 28 ) 

k k 
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where ojq is the energy difference between the excited 
state |1) and the ground state |0) of the qubit. Bosonic 
creation and annihilation operators bj, and 6 ^ correspond 
to the fcth bath mode with frequency cuk, and gk are the 
coupling constants. 

Suppose that at time t = 0 the state of the composite 
system (S+B) is characterized by some density matrix 
Pg^(O). Then at time t the average value of a Heisenberg 
picture operator A(t) is given by 

(A(t)) =Tr{exp(ifft)Aexp(—iHt)gg^(0)}. (29) 

Below, the notation (H) will be used for averages at f = 0. 

The quantities of principal interest are the coherences 
(CTj_(t)), where = {ai ±*< 72 ) /2. They are related di¬ 
rectly to the off-diagonal elements of the reduced density 
matrix of the qubit: 

{cr+{t)) = (0|P5(<)|1), (cr_(t)) = (llesWIO). (30) 

where 

SsW = Tr^ {exp{-iHt)gsgiO) exp{iHt)} . (31) 

The dephasing model (1^ has two distinctive features. 
First, the operator Ug commutes with the Hamiltonian 
and, consequently, the average populations of the canon¬ 
ical states (HB do not depend on time. Thus we have a 
unique situation where the system relaxation may be in¬ 
terpreted physically as “pure” decoherence and exchange 
of entropy [o, [3 rather than dissipation of energy. Sec¬ 
ond, in this model the equations of motion for all releyant 
operators can be solyed exactly 0. This allows one to 
study the time eyolution of the coherences for different 
initial conditions. Here we leaye out many details for 
which we refer to Ref. and simply quote some impor¬ 
tant results. 


If the initial state is prepared by a non-selectiye mea¬ 
surement, then, taking the initial density matrix of the 
composite system in the form o, we get for the coher¬ 
ences (1501) 

{cr±{t)) = [H)„cr±(f)H„e"^^] , (32) 

m 

where /3 = l/fcgT, and Z = Tr {exp(—/3H)} is the equi¬ 
librium partition function. The analogous formula for the 
case of a selective measurement [see Eq. ©] is evident. 

As shown in Ref. Q, the time-dependent qubit opera¬ 
tors crj_ (<) in the dephasing model (1051) are given by 

a± (t) = exp {±iuJot =F R{t)} a± (33) 

with 

m = E 

k 

(34) 

Using the above expressions and the exact relations 

g-/3ff|0) = g/5‘^o/2g-^ffB ' (g) |Q\ 

(35) 

g-/3ir|i) _ g-/3i^o/2g-/3ffB ' (g) |1)^ 

where 

H^b'^ =Y,^khlbu±Y.^9khl + glh^). (36) 

k k 

it is a straightforward matter to carry out the trace over 
the bath degrees of freedom in Eq. (l32l) . After some 
algebra (for details see Ref. i) , one obtains 


\akit)bl -al{t)bk 


ak{t)= 2gk 


X _ 
W/c 


J 


{cr±it)) = (cr±)e 


J2m |(0|HLcr±H„i|0)e(^“o/2 -I- (l|HLcr±H„j|l)e“(^“o/2| 


(37) 


with the initial coherences 


of Eq. ((57)) . we have 


(c^±) 


1 

2 cosh (/3 wq/2) 




+ (38) 


In the case of a selective measurement when the qubit 
is prepared in a pure quantum state {ip) = CglO) -I- Cgjl) 
with |cqP -I- jcgp = I, the initial density matrix of the 
composite system is taken in the form ([0|). Then, instead 


(a±(t)) = (a±)e±““‘e-^(*) 

|gp|2g/3wo/2±i<I.(t) _|_ |^^|2g-/3a;o/2=FJ$(t) 


(39) 


where (crg_) = {tplcj-lip). 

Formulas (157)) and (l39l) contain two relevant functions. 
The so-called decoherence function 7 (f) is defined as 


7 (f) = / dw J(a;) coth(/3w/2)-, (40) 

Jo ^ 










5 


where the continuum limit of the bath modes is per¬ 
formed, and the spectral density J(w) is introduced by 
the rule 


_ poo 

'^M9k\^fi^k)= dujJ{uj)f{uj). (41) 

k -^0 

It is clear that 7 (<) is precisely the quantity which de¬ 
termines the relaxation of the off-diagonals (l3^ due to 
vacuum and thermal fluctuations in the bath [l0[. The 
other function, 4>(t), is given by 

$(t)= rdujJiLof-^. (42) 

Jo ^ 

As discussed in Ref. [^, this function accounts for the 
influence of initial qubit-environment correlations on the 
dynamics of decoherence. These correlations are inher¬ 
ited from the pre-measurement equilibrium state due to 
the presence of the interaction term in the total Hamil¬ 
tonian (E51) . Mathematically, it has the effect that the 
operators Hg and Hg [see Eqs. (l35|) and (l36l) ] differ in 
the sign of the interaction term. 


is the effective decoherence function including the corre¬ 
lation contribution 


r(j)= “2 In 


1 - 


(l-(a3)2)sin2$(t) 


[cosh(/3wo/2) - (cr 3 )sinh(/ 3 wo/ 2 )] 2 j ’ 

(45) 

whereas x(J) is the time-dependent phase shift with 


tanx(t) 


sinh(/3wo/2) - ( 0 - 3 ) cosh(/3wo/2) 
cosh(/3wo/2) - {as) sinh(/3u;o/2) "" ^ ’ 


(46) 

In writing the above formulas we have used the obvious 
relations ( 0 - 3 ) = {ip\as\Jj) = [cip - jcQp. 


It is important to note and easy to see from Eqs. go]) 
and (HSl) that both terms in the effective decoherence 
function are always positive. Thus, in cases where 
the initial state is prepared by a selective measurement, 
initial qubit-bath correlations may be viewed as an addi¬ 
tional source of decoherence. 

C. Non-selective measurements 


B. Selective measurements 

The result (15^ for selective measurements has been 
analyzed in detail in Ref. Q, so that here we merely 
touch on some relevant points. 

The expression (1391) can be rewritten more transpar¬ 
ently as 

{a±{t)) = (cr±) exp[±i(a;ot -f x(t)] exp[- 7 (t)], (43) 

where 


Now we turn to Eq. dSZl) and consider the general non- 
selective measurement scheme described by Eqs. m- 
Using the representation m, we obtain after some alge¬ 
bra the expression (1431) in which the correlation part of 
the decoherence function is now given by 

lcor{t) = “2 jj2 - ij sm2$(t) 

+ §sin(2$(t))|, (47) 


lit) = lit) + Icorit) 


(44) where we have introduced 


Ni = sin^(6(a/2) - cos^(6(a/2)} sin^ 6»i -h cos^{9a/2) - sin^(6»a/2)} sin^ 62 

+ sinh(/3wo) sin^ 6 a cos sin 9i sin 02, 

^ 2 = 2 cos0a sinA^ sin0i sin02, (48) 

74 = {4 sin^ 9a + sm\9a/2) + cos4(0a/2)} sin^ 0i + {i sin^ 9a + cos4(0a/2) + sin4(0a/2)} sin^ 02 

-I- {cosh(/Iwo) sin^ 0a + 2 [sin"‘(0a/2) -|- cos'^(0o/2)] } cos A^ sin0i sin 02 


To simplify notation, we have written 0^, (|)^ for 0^., 0^., 
and denoted A,^ = 4>i~ 4’’^.- T4ie expression for the phase 
shift xit) is 


x(t) = arctan 


Ni sin$(t) 


Dcos^{t) + N 2 sin$(t) 


(49) 


The initial coherences (l38l) can be directly evaluated to 
yield 


(o'±) = 




4cosh(/3a;o/2) 

X •|sin0j^ 

-(-gT*A,^ sin 02 


2 2 


cos ~ 2 ~^^ 


sin 


.(50) 
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Formulas (gzl) - dMl), together with Eqs. dMl) and (l4^ . 
determine the time-dependent coherences (ctj_(<)) or, 
what is the same, the off-diagonals of the qubit density 
matrix [see Eq. (15(11) ]. It is often convenient to use the 
representation of the qubit densit y rn atrix gg(t) in terms 
of the Bloch vector v{t) = {(T{t)) 

esW = + (SI) 

The magnitude of the Bloch vector satisfies 0 < u(t) < 1, 
and V = 1 only if the qubit is in a pure quantum state. 
It is easy to check that 

V{t)= [4{cr+(t))(cr_(t)) -h (52) 

Since commutes with the Hamiltonian (j28p . we have 
(tT 3 (t)) = (cts). Taking also into account Eq. (|43)) . we 
obtain 


u(t) 


4{(7-|_)(cr_)e 


+ 



(53) 


Thus, to calculate v{t), all we need is the initial average 
((T 3 ). For the general non-selective measurement scheme 
m, this average can be found by using the formula anal¬ 
ogous to ((551) . A straightforward algebra gives 


(<73) = 


2 cosh(/3wo/2) 
X { cos 9i 






+ cos 6*2 


e/3-o/2^0g2^^g-/3.;o/2gi^2^ 


■ (54) 


D. Some special non-selective preparation 
measurements 


To gain an insight into new features of the qubit dy¬ 
namics in cases when the initial state is prepared by a 
non-selective measurement, we will apply the above gen¬ 
eral expressions to some special preparation schemes de¬ 
scribed in Sec. im 

We start with the scheme (1251) for which, according to 
relations m, 


6 * 1 -I- 6*2 = TT, sinA^ = 0, cosA^ = —1, (55) 


where 9^ = 0^ and (j)-^ = Noting that in this case 
^2 = 0, Eq. (l47l) is manipulated to the simple form 


'Jcorii') — 2 


1 -H 


sin^ $(t) 


sinh"^(/3a;o/2) 


For the phase shift (15^ we find 

tanx(t) = coth(/3a;o/2) tan$( 6 ). 


(56) 


(57) 


The results (1551) and (1571) have several notable proper¬ 
ties. First, they are universal in the sense that they do 


not depend on the qubit states ja) and |6) in Eqs. (l23l) 
describing this type of non-selective measurements. In 
particular, the same expressions for jcorit) and x(t) hold 
for the simplest scheme ((22l) where \b) = ja), i.e., the 
measuring device does not disturb the basis states. Next, 
we note that the function (1551) satisfies ^cor{t) < 0 at all 
times t. In other words, we have an enhancement of 
coherence in the qubit caused by initial qubit-bath cor¬ 
relations! Moreover, it is seen from Eq. (1551) . that the 
l7cor(6)| grows with temperature, so that in the temper¬ 
ature range /Swg <C 1 the effective decoherence function 
(155)) may even become negative, at least during the initial 
stage of the system’s evolution. To illustrate this point, 
we have evaluated the reduced coherence 

Ko'(0)l/l(c^)l = l(CT±(t))|/|(cr±)| = exp [- 7 ( 6 )] (58) 

using Eqs. (|55)). (|55)). and dSH). The bath spectral density 
was taken in the form 

J(a;) = A« (59) 

which is most commonly used in the theory of spin-boson 
systems [Sill 13 HI. Here lvc stands for some “cutoff” 
frequency, and As is a dimensionless coupling constant. 
The “dynamical part ” (gOj) of the decoherence function 
and the function 4)(6), Eq. (1551) . have been studied in 
detail in Ref. Q for the sub-Ohmic (0 < s < 1), Ohmic 
(s = 1), and super-Ohmic (s > 1) cases. Here we shall re¬ 
strict ourselves to the most prominent Ohmic case where 

4>(6) = A arctan(a;c6), A = Ai, (60) 

and the time behavior of the correlation term (1561) is very 
sensitive to the value of the coupling constant A. 

Figure [T] shows the time dependence of the reduced 
coherence in the temperature range /3u}q <C 1. At the 
initial stage of evolution {ojJ: < 1) the correlation effects 
dominate, so that the decoherence function dSI) takes 
negative values. At times ujJ, ^ 1 the correlation effects 
are suppressed due to vacuum and thermal fluctuations 
contributing to the “dynamical part ” (gni) of the deco¬ 
herence function. Using the explicit expressions for the 
function ((551) in the super-Ohmic (s > 1) and sub-Ohmic 
(0 < s < 1) cases (see Ref. Q), it can be shown that 
a similar time behavior of the reduced coherence is ex¬ 
pected in these coupling regimes. 

It is worthwhile remarking that even for a moderate 
qubit-bath coupling (A ~ 1), the maximum value of the 
coherences |((73_(6))| may be by one order of magnitude 
larger than the initial value |((T 3 _)|. This fact appears at 
first rather strange and even paradoxical. In particular, 
it seems likely that for some initial conditions, the mag¬ 
nitude of the Bloch vector (1551) may exceed the unity. 
This is, of course, not so for the following reason. The 
point is that the initial averages (1551) and ((55() depend 
on temperature. For the preparation procedure under 
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FIG. 1: Time evolution of the reduced coherence in the 
Ohmic case for different values of the dimensionless cou¬ 
pling constant. The preparation measurement is described 
by Eqs. (I23II . Parameter values: Pujq — 0.1, uJojoJc = 0.01. 


consideration [see Eqs. m], we find 

(o-i) =-^e=^*‘^'>tanh(/3wo/2)cos6»„sin6»^, 

(61) 

(o's) = — tanh {Pujo/2) cos 9^ cos 0 ^. 

In the temperature range ^ujq <C 1, we have |(crj_)| ^ 
1 and KcTa)! <?C 1. Recalling Eqs. (1531) and (1^ , one 
can show that, at all times, v{t) < 1 , as it must be. 
Other physical questions related to the enhancement of 
coherence for this preparation scheme will be discussed 
in Sec. ITVl 

Now we shall consider the dynamics of decoherence 
in the case of another notable non-selective preparation 
scheme described by Eqs. (EH). This scheme corresponds 
to the following choice in Eqs. (H51) : 


6 »i = 6»2 = 6 »f„ sinA. = 0, cosA. = 1 . (62) 


Erom the second of Eqs. (l48l) it is clear that N 2 is again 
zero. Then, after some simple algebra Eqs. 611) and 63) 
give 


Icorii) = - - In 


1 - 


sin^ $(t) 


cosn (/3wo/2) 


(63) 


tanx(t) = tanh(/3a;o/2) tan$(t). 


(64) 


The initial averages (1501) and ((50)) now take the form 


(o'i) = sin 6 (j^, (cts) = cos9f^. (65) 


Similar to formulas (1501) and (1571) . the results (1001) and 
( 1001 ) are universal in the sense that they do not depend 
on the qubit states |a) and \b) in Eqs. (1051) . Another 
distinctive property of Eq. OMl) is that 7 cot -(0 > 0 at all 


times. Physically, in this case the initial qubit-bath corre¬ 
lations lead to additional decoherence. It is also interest¬ 
ing to note that expressions ( 1001 ) and ( 1001 ) are identical to 
expressions 63 and 63 for the selective measurement 
with ((Tg) = 0 , i.e., with equal populations of the basis 
states 63- This is not accidental; for a discussion see 
Appendix [^ 

A comparison of Eqs. (1551) and (1021) . together with 
the fact that the Euler angles 9^ and ^2 enter the func¬ 
tions (1051) only through sin 0 g and sin 02 , suggests that 
A^ = '/'i ~ </'2 is a key quantity determining the main 
qualitative features of the system’s evolution. To illus¬ 
trate this point, suppose that the system is initially pre¬ 
pared by the non-selective measurement (O with the 
Euler angles of the qubit states | 6 g) and [ 62 ) satisfying 

6 »i = 02 , sinA^ = 0, cosA^ =-I. ( 66 ) 

The corresponding initial averages (1501) and (1501) are 

(cTj.) = — tanh(/3a;o/2)cos0Qsin6(i, 

(ct3)=cos0i. 

We mention that formulas (1001) differ from Eqs. (1001) only 
in that cos A^ is now of opposite sign. Nevertheless, it is 
evident that we obtain for 7 cot-( 0 the result (1501) which 
corresponds to the entirely different evolution of the co¬ 
herences ((Tj_(t)). 

It would be instructive to look at the last example 
from another point of view. Let us write the post¬ 
measurement qubit state |5g) in the canonical basis dill) : 

1 ^ 1 ) = co|0)-bci|l), (68) 

where the amplitudes Cg and Ci can be expressed in 
terms of the Euler angles by using Eq. 63- Then us¬ 
ing Eqs. dMl) leads to the following representation for 
the state I 62 ): 


I 62 ) =i(co|0)-ci|l)). (69) 

There is no new physics in the appearance of i, but the 
additional phase shift between the basis states, as com¬ 
pared to Eq. (1551) . radically influences the qubit’s dynam¬ 
ics. 

Thus far we have been concerned with special types of 
non-selective measurement schemes which lead to physi¬ 
cally important features of decoherence. To give a com¬ 
prehensive review of all possible regimes of evolution, one 
should appeal to Eq. (HTl) . In general, the quantities (ITTl) - 
(l49l) are rather complicated functions of the polar angles 
9g^, 9-y, 02 , and the difference A,^ = (j>-^ — (j) 2 of the az¬ 
imuthal angles. In addition, they depend on temperature 
and the qubit energy Wg. This makes a detailed analysis 
of 7cor(0 rather cumbersome for our discussion. Nev¬ 
ertheless, we can formulate a simple sufficient condition 
that the initial preparation of the system by a general 
non-selective measurement pO) leads to enhancement of 



















coherence in the qubit. Notice that for all qubit states 
|5^) and I& 2 ) with sinA^ = 0 , we have N 2 = 0 and, con¬ 
sequently, the term with sin(2$(f)) in Eq. (lETll vanishes. 
Then we arrive at the conclusion that Jcori't) < 0 at all 
times if 

sinA,^=0. (70) 

Clearly, the second condition implies cosA^ = ±1. We 
will not give here a somewhat lengthy formal analysis 
of Eqs. (1701) since it does not add anything substantially 
new to the results of the above discussion. We only note 
that the enhancement of coherence takes place for the 
post-measurement states |&]^) and I 62 ) with cosA^ = — 1 
and + 02 ss tt or 6*2 ~ ^ 2 - other words, the post¬ 
measurement states should be close to the states in the 
measurement schemes HMD or (1551) . 

IV. THE PURITY AND ENTROPY OF THE 
QUBIT 

If at time t the effective decoherence function is 
negative, the state of the qubit is “less mixed” than ini¬ 
tially. This property can be characterized quantitatively 
by the von-Neumann-Shannon information entropy 

S{t) = -Trg{£>g(f) Inggit)}. (71) 

Using Eq. the information entropy of a qubit can 

be expressed in terms of the Bloch vector magnitude Q : 

S{t) = ln2 — i (1 -I- u) In (I -|- u) — i (1 — w) In (1 — u). 

(72) 

Since 0 < u < 1, we have 0 < S' < In 2 with S = 0 for a 
pure quantum state (u = I). 

Another measure of the “mixedness” (or the lack of 
information about a system) is the so-called purity of 
the system’s state : 

Vit) = Tys{gUt)}. (73) 

Again using Eq. (1511) . we obtain for a qubit 

nt) = ^{l + v^). (74) 

Obviously 1/2 < P < I with P = 1 for a pure state. 

In principle, either S{t) or V{t) may be used to mea¬ 
sure the degree of coherence in a qubit. In both cases the 
key quantity is the magnitude of the Bloch vector ([53|) . 
Here we shall discuss the time behavior of the purity and 
entropy in the non-selective measurement schemes (1551) 
and (I 66 p for which the enhancement of coherence is most 
pronounced. 

Figures [5] and [3] illustrate the evolution of the qubit 
entropy d72l) and purity d74l) at a fixed temperature for 
different values of the coupling constant. 

It is seen that with increasing the coupling constant 
the maximum purity becomes larger and is shifted to 



FIG. 2: Time evolution of the qubit entropy S{t) in the Ohmic 
case for different values of the coupling constant: A = 2 (B 
and E), A = 4 (C and F), A = 6 (D and G). Filled sym¬ 
bols correspond to the preparation measurement described 
by Eqs. dSSJ, half-filled - to the measurement described by 
Eqs. (1661) . Other parameter values: Puuo = 1, Loojujc = 0.1, 
Oa = 0, 6 I 1 = 7r/4. 



FIG. 3: Time evolution of the qubit purity P(t) in the Ohmic 
case for different values of the coupling constant. The symbols 
and the system parameters are the same as in Fig. [T] 


smaller t. For large values of the coupling constant, there 
appear oscillations in both the purity and entropy (cf. 
Fig. [T]). According to Eqs. (I5T]) and (l67l) . the initial 
coherences |(cr£)| are identical in the preparation mea¬ 
surement schemes described by Eqs. (1551) and (1551) . Note 
also that the time dependence of the decoherence func¬ 
tion 7 (t) is the same in both cases. The corresponding 
pairs of lines in Fig. [3] are therefore just shifted vertically 
from each other because of the different values of (< 73 ). 
Expressions dH]) and (1571) for {a 3 ) show that, at all tern- 
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peratures, the preparation scheme (l5^ leads to less pure 
states of the qubit as compared to the scheme (1551) . 

Figures Eland [5] display the evolution of the qubit pu¬ 
rity and entropy at different temperatures but for a fixed 
value of the coupling constant. 



FIG. 4: Time evolution of the qubit entropy S{t) in the Ohmic 
case at different temperatures: fitoo ~ 0.01 (B and E), Puuo = 
0.1 (C and F), Pljq = 1 (D and G). Filled symbols correspond 
to the preparation measurement described by Eqs. (1551) . half- 
filled - to the measurement described by Eqs. (1661) . Other 
parameter values: A = 6, wo/cUc = 0.1, 6 a = 0, = 7r/4. 



FIG. 5: Time evolution of the qubit purity P(t) in the Ohmic 
case at different temperatures. The symbols and the system 
parameters are the same as in Fig. (4] 

It is seen that lowering of temperature is favourable 
for purification of the qubit state. To explain this, it 
should be remembered that the time evolution of the pu¬ 
rity V{t) and the entropy S{t) is the net result of two 


processes. First, there is the enhancement of coherence 
due to the combined effect of the non-selective prepara¬ 
tion measurement and qubit-bath correlations. Second, 
thermal fluctuations in the bath suppress the coherence 
effects. At high temperatures, the latter process domi¬ 
nates. Note in this connection that the maximum of the 
purity shifts to smaller t with increasing temperature (see 
Fig.[S]). 


V. CONCLUSIONS 

In this paper we have examined the properties of the 
reduced qubit dynamics in cases where the initial state 
of the composite system (qubit plus environment) is 
prepared through non-selective quantum measurements. 
Our main result is that for some preparation schemes 
the interplay of the measurement process and qubit- 
environment correlations can lead to a significant en¬ 
hancement of coherence in the qubit during the initial 
stage of evolution. The non-trivial feature of this effect 
is that, in general, a non-selective measurement produces 
a mixed state of the composite system. Then the purity 
of the qubit’s state grows while its entropy decreases with 
time until thermal fluctuations suppress this environmen¬ 
tally induced “purification” process. It deserves to be 
pointed out that the temperature dependence of the pu¬ 
rity growth is determined by several factors. First, the 
initial coherences (i.e., the off-diagonal elements of the 
qubit density matrix) generated by a non-selective mea¬ 
surement decrease dramatically with increasing temper¬ 
ature. Second, the destructive effect of thermal fluctu¬ 
ations in the bath also becomes important just at high 
temperatures. But surprisingly, the dynamical enhance¬ 
ment of coherence due to the qubit-bath correlations 
grows with temperature, so that the maximum of the 
coherences |(crj_(t))| may be much larger than the initial 
values |(cr±)| (see Sec. IIII|) . Due to the interplay of the 
above factors, the resulting purity of the qubit states de¬ 
creases with temperature, but not so rapidly as one might 
expect from intuitive considerations. 

Summarizing, if the initial state is prepared by a non- 
selective measurement, it is possible to achieve the en¬ 
vironmentally induced purity growth in a qubit during 
some time interval by setting the measurement device in 
a proper way. This effect may be of interest in view of its 
connection with problems of measurement-based quan¬ 
tum control in open quantum systems. We refer, e.g., to 
the recent paper were the dephasing model (l28l) was 
used to study the influence of system-environment corre¬ 
lations on the so-called quantum Zeno and anti-Zeno ef¬ 
fects in repeated selective measurements on single-qubit 
and many-qubit systems. 

As a final remark we wish to emphasize that our study 
of the environmentally induced purity growth does not 
claim to be complete even on the single-qubit level. Here 
we will touch briefly on two open questions which de¬ 
serve further investigation. First, our analysis was based 
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on the assumption that initial correlations between the 
qubit and the environment are inherited from the pre¬ 
measurement equilibrium state due to the presence of 
the interaction term in the total Hamiltonian H. This as¬ 
sumption is adequate in describing many real situations 
and is commonly accepted in the theory of open quan¬ 
tum systems lH-il. Nevertheless, one may imagine a 
variety of different correlated initial states by replacing 
Pgq in Eq. (HD by some nonequilibrium premeasurement 
density matrix g'. One should note, however, that within 
such a general formulation of the problem the trace over 
the bath degrees of freedom in {a^ {t)) cannot be carried 
out without a detailed information about physically rea¬ 
sonable forms of g'. The second important point is that 
the model ((28)l describes only the dephasing mechanism 
of decoherence. Although this mechanism can dominate 
in real physical systems [13, HI] , in general the popula¬ 
tion decay should be taken into account. However, it is 
a challenging problem because in this case the model is 
no longer exactly solvable. 


Appendix A: Preparation of pure qubit states 
through non-selective measurements 

Here we briefly discuss the connection between the 
non-selective preparation scheme (1^ [see also Eqs. (15^ ] 
and selective measurements. 

Using the explicit expressions for the H-operators, the 
initial density matrix (HD of the composite system is writ¬ 
ten as 


esBiO) = \b){b\<^QB{a) (B.l) 

with the bath density matrix 

gB{a) = (d|£<eq|a) + (““beql “ «) ■ (B.2) 

Eormula (IB. II) shows that after a non-selective measure¬ 
ment of this type, the qubit is prepared in a pure state 
16). Note also that the qubit and the bath are completely 
uncoupled since the bath state does not depend on the 
qubit state, and vice versa. Although the bath density 
matrix dEl formally depends on the basis state ja) de¬ 
termining the effects in Eqs. (EH), it is easy to see 
that ggia) in fact is independent of ja). Indeed, since the 
states |d) and j — a) form an orthonormal basis, Eq. (IB.2|) 


may be rewritten as £* 3 ( 0 ) = Trsp^^. If so, the trace can 
now be calculated with any other orthonormal basis. In 
particular, it can be done with the canonical states HID, 
so that the initial bath density matrix takes a universal 
form 


Qb = (O|0eq|O) + (ll^eqll)- (B.3) 


Let us now assume that the qubit is initially prepared 
in some pure state {ip) through a selective measurement. 
Then the initial state of the composite system is given 
by formula (HD or, what is the same, by 


0Sb(.O) = ^ q'bW 

with the initial density matrix of the bath 


QbW 


(V^lgeqlV’) 

TrB(V^leeqlV’)' 


(B.4) 


(B.5) 


In this case the bath carries information on the qubit 
state since its density matrix depends on \ip) through the 
interaction term in the Hamiltonian. Consequently, the 
product (IB. 41) should be interpreted as a correlated state 
of the composite system Q. For the dephasing model 
(I28L the bath density matrix (IBAl) can be written in 
a more transparent form. To do this, we note that the 
equilibrium density matrix g^^ of the composite system is 
diagonal with respect to the canonical qubit states HD- 
Then, writing lip) as a decomposition \ip) = Cq|0) -|-C;^|1), 
we obtain for the density matrix (IB31) : 


|CoP(0|£^eq|0) + |CiP(I|Peq|l) 

|coPTr^(O|0eq|O) + |ciPTr5(l|£.eq|l)' ^ ^ 

Since, in general, the bath density matrices (jED and 
(IBAI) differ from each other, the time evolution of the 
composite system depends on the type of the prepara¬ 
tion measurement. Suppose, however, that the qubit is 
prepared selectively in a pure state \ip) with equal pop¬ 
ulations of the canonical basis states, i.e., with (a^) = 
{ipla^lpj) = 0. In this case we have |cqP = \ci\^, so that 
the density matrices (IB.31) and (IB.61) coincide. This is 
the reason why the decoherence dynamics for the non- 
selective measurement scheme (1^51) is identical to the 
decoherence dynamics for a selective measurement with 

(cTa) = 0. 
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